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ABSTRACT 


Synthetic aperture sonar beamforming and signal proc- 


essing relies on properly steering and focusing the aper- 





ture beam pattern in order to co-phase all the received 





Signals. Due to the effects of motion in the synthetic 





aperture sonar problem, the propagation path between the 


transmitter, discrete point scatterer, and the receiver is 





time varying. Traditionally, simple approximations are 











used to determine these propagation ranges and angles of 











incidence and scatter. Methods to determine these ranges 


and angles exactly may significantly improve array gain 








and, therefore, target detection. 


This thesis investigates improvements to SAS signal 





processing algorithms using exact methods for the calcula- 
tion of the time-varying ranges between transmitter and 


discrete point scatter, and between discrete point scatter 





and receiver, and the phase angle of the scattered acoustic 
signal incident upon the receiver. Using computer simula- 


tions, exact range and angle calculations were performed 





for different scenarios and compared to ranges and angles 


determined using standard approximations. The exact ranges 








were then used to determine incident phase, and were again 
compared to the approximate methods. Comparison of the 


exact and approximate methods was based on range estimation 














error and percentag rror. Improvements in synthetic 


aperture array gain using exact phase weights based on 








exact, time-varying range solutions are proposed. 
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EXECUTIVE SUMMARY 


The goal of this thesis was to demonstrate improve- 
ments to Synthetic Aperture Sonar (SAS) beamforming algo- 


rithms. SAS systems are currently in development to sup- 





port the Navy’s Long-Term Mine Reconnaissance System 
(LMRS). The LMRS will use a submarine-launched autonomous 


underwater vehicle to create a synthetic array for high- 





resolution mine imaging and classification. 











Currently, SAS signal processing techniques make vari- 


ous assumptions and simplifications in the determination of 








the signal propagation ranges and phase shift corrections 

required for accurate acoustic modeling. These approximate 
techniques simplify the steering of the beam pattern of the 
sonar aperture but also unnecessarily introduce errors that 


affect image quality. 





Beamforming techniques have been developed that can 








exactly determine the parameters required to steer the beam 

















pattern of a sonar array. Using this exact method, acous-— 





tic signal propagation in the ocean for a SAS system was 








simulated. This simulation was conducted using stationary 











targets to replicate the problem of water-borne mine detec 








tion and classification. Water-borne mines include buoy- 





ant, tethered, and bottom ocean mines. The primary focus 


of these simulations was to demonstrate the improvement in 





SAS beamforming that can be realized through exact propaga- 


tion range calculation. With exact range estimation, exact 





phase weights are determined and used to improve array 


gain. 


x1il 


An analysis of the limitations of the most common ap- 


proximations is presented. The exact beamforming methods 





have no restrictions on their application, and this is 


demonstrated through simulation. The simulation of the 





propagation of signals from transmitter to target and then 





to the receiver allows direct comparison of the standard 





approximate methods with the new exact methods. Estimation 


of signal propagation range, travel time, angles of inci- 





dence and scatter at the target and angles of incidence at 








the receiver are all improved. Finally, fast Fourier 


transform beamforming is conducted and the performance of 





the exact and approximate methods is compared using a sim- 





ple signal composed of three different sinusoids. 





This research demonstrates that simple approximation 








techniques traditionally used to estimate SAS propagation 


ranges are not useful in current SAS projects that employ 





high frequency signals and long target ranges. This re- 


search also shows that SAS propagation ranges can be accu- 





rately approximated using the binomial approximation propa- 








gation model presented in this thesis. Although the bino- 





mial approximation propagation model is accurate under most 





conditions, this model places some restrictions on target 


range and platform speed. The effects of operating outside 





the binomial approximation model limits are demonstrated. 











It is shown that implementation of the exact propagation 





range method does not significantly add to the complexity 

















of the beamforming algorithm. The exact method introduces 








no error and places no restrictions on target location and 





platform speed. Use of the exact propagation model is 
recommended to allow maximum accuracy and flexibility as 


SAS systems evolve. Further investigation will determine 





Xiv 


the improvements that can be realized by applying the exact 


time-varying angles of incidence and scatter to the evalua- 








tion of target scattering functions. 
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na INTRODUCTION 


This chapter provides an overview of synthetic aper- 
ture sonar concepts. The modeling framework used to simu- 


late signal propagation through the ocean medium is also 





introduced. 


A. SYNTHETIC APERTURE SONAR 
Synthetic Aperture Sonar (SAS) systems attempt to ex- 


ploit the benefits that can be obtained from very long 





sonar arrays while maintaining the actual physical size of 





the array to be relatively small. This is an attractive 





technique because it would allow deployment of Unmanned 








Underwater Vehicles (UUV) to detect, localize and classify 


mines in the littoral environment. Vehic] size is limited 











by the launching capability of the controlling ship or 


aircraft and a need to conduct covert surveillance and mine 





hunting. The high resolution necessary for classification 





of mines, however, requires a very long array. One of the 





most attractive ways to satisfy these competing require- 








ments is through aperture synthesis. 


The principle of aperture synthesis is the coherent 


combination of successive returns from a transmitter lo- 








cated on a moving platform. This coherent processing of 








received echoes provides enhanced resolution in azimuth 


compared with standard sonar systems. The key to SAS sig- 





nal processing is estimating the time delay of each re- 
ceived pulse while compensating for motion of the array and 


the target. This process can be thought of as focusing of 


the phase history of the signal [1]. Traditional sonar 





arrays use multiple transmit and receive elements to simul- 
taneously transmit and then simultaneously receive acoustic 
Signals. Conversely, synthetic aperture arrays synthesize 


the array by processing the output electrical signals of a 





single transmit and receive element or group of elements as 





the element travels through the medium. A comparison of 
the traditional array and a synthetic aperture can be seen 


in Figure 1. 









PING 1 






simultaneous 
transmission 
& reception 


sequential 
transmission 
& reception 


Traditional Array Synthetic Aperture Array 


Figure 1. Traditional and synthetic aperture arrays 
[After Ref. 2]. 


SAS systems are similar to Synthetic Aperture Radar 
(SAR) systems in many respects but SAS signal processing 
provides some unique challenges. SAR systems have wider 
absolute bandwidth and higher carrier-to-bandwidth ratio, 
or quality factor QO. The high Q nature of SAR dictates 
that typical motion and range estimation errors have little 


effect on the echo. The errors result ina shift of a few 





centimeters on a signal envelope of a few meters. In con- 


trast, SAS motion and range estimation errors of a few 








centimeters result in a shift of the target envelope that 


2 





is also on the order of a few centimeters [3]. The strong 
dependence of SAS systems on accurate range and motion 
estimation is the focus of the research presented in this 


thesis. 





Several other challenges must be considered when d 
signing a SAS system. The speed of the platform must be 


high enough to maintain control of the submerged vehicl 





but slow enough to prevent undersampling of the target 





area. As described above, SAS system resolution can be 


significantly degraded by residual motion errors of the 





transmit and receive platforms. 





B. THE OCEAN AS A LINEAR, TIME-VARIANT, SPACE-VARIANT 
FILTER 


Computer simulation of acoustic ray propagation in an 





unbounded ocean medium was the primary means of evaluating 
the exact propagation equations and comparing these tech- 


niques to common approximations. The mathematical frame- 





work used for these simulations models the ocean as a lin- 
ear, time-variant, space-variant filter. The theory, meth- 
ods, and derivation of this linear systems approach to 


acoustic field theory are contained in [4]. A block dia- 





gram model of pulse propagation is shown in Figure 2. The 


Simulations conducted in this thesis focus on determination 








of exact and approximate acoustic propagation paths in- 





cluded in the complex frequency response of the ocean m 
dium. The complex frequency response of the ocean medium 


at time ¢, position r, and frequency f due to application 


th) 1 as 








Of “a unit-amplitude impulse at position 7 is H,, (tr 


seen in Figure 2. 














xytn) | Ayltrl4.m) | y,@r) 


H,dt, rlf, r;) ¥,,03,9r) 


xr) 





Af.) 






X(f,0) Yr) 





Figure 2. Block diagram model of small-amplitude pulse 
propagation [From Ref. 5]. 


Models used for simulation are based on bistatic scat- 


tering with a stationary target. This scenario is intended 











to model the performance of a submerged vehicle using a 


side scan sonar to image stationary mines. The transmitter 








and receiver are co-located on the platform. Due to motion 








of the vehicle the position of the transmitter and receiver 

















is time varying. It is important to note that the vehicle 


will be in motion during the entire transmit and receive 








cycle. Some SAS signal-processing algorithms assume that 
the vehicle is stationary during each transmit and receive 
cycle and moves only between each cycle [6]. This simpli- 
fication requires motion compensation to be applied during 
signal processing and can introduce errors in the image. 
The models used in this thesis are exact and require no 


motion compensation. The exact models do, however, assume 





that the velocity of the platform is constant during trans-— 





mission and reception of individual pulses. 


Ci RESEARCH GOALS 
This thesis focuses on the direct comparison of exact 
and approximate signal propagation models as applied to a 


SAS system employed to image stationary targets. This 








research validates, through computer simulation, the exact 





propagation model. The exact model is then used to evalu- 





4 


ate errors introduced by approximate methods. The goal of 


this research is to provide a proven SAS signal-processing 





algorithm that removes errors due to propagation range, 
angle, and platform motion estimation. With reduced error, 


the synthetic aperture array has increased array gain and 





improved image quality. 


D. DOD RELEVANCE 


Mine detection and classification is one of the most 








important challenges faced by the Navy. Simple, inexpen 





sive, effective mines are available to all our potential 








adversaries. Deployment of water-borne mines by an adver- 





sary precludes the introduction of ships or troops into an 





otherwise unprotected area. 





Th ffectiveness of mine warfare was proven in the 


1991 Persian Gulf War. The USS Tripoli sustained 3.5 mil- 





lion dollars in damage after being struck by a 1500 dollar 


mine. Later in the war, the USS Princeton was disabled by 





a Similar mine. Repairs to the USS Princeton cost 24 mil- 


lion dollars [7]. Battle group commanders are justified in 











their reluctance to employ assets where the threat of wa- 








ter-borne mines exists. 





To address the asymmetric threat of mines, the U.S. 
Navy is currently developing a Long-Term Mine Reconnais-— 
sance System (LMRS). The LMRS is an unmanned underwater 
vehicle that will use a forward-looking sonar system for 
mine detection and a side-scan sonar system for localiza-— 


tion and classification. The side-scan sonar is used as 











the transmitter and receiver for a SAS system. This system 


provides no risk access to mined waters while collecting 





information essential to mine hunting operations. A con- 


ceptual drawing of the LMRS can be seen in Figure 3. 





Figure 3. LMRS conceptual drawing [From Ref. 8]. 


E. SUMMARY 


Reliable detection and classification of water-borne 








mines is an essential capability for the U.S. Navy. SAS 








systems on UUV’s are currently under development to address 
this challenge. Traditional SAS beamforming algorithms use 
approximations to calculate acoustic signal propagation 


ranges and angles of incidence. Exact methods of determin-— 





ing these parameters have been developed and improve system 


performance. This thesis uses computer simulations of 





exact and approximate propagation models based on linear 


systems theory to compare the performance of these methods. 


F. THESIS ORGANIZATION 


This thesis is organized into five chapters. Chapter 





provides an overview of the SAS concept, military rele- 





vance of this research, and a general description of the 


modeling methods and evaluation techniques used for system 





simulation. An overview of the application of linear sys- 


tems theory to acoustic propagation modeling is presented 


























in Chapter . Chapter introduces the mathematical 








models used for simulation of SAS as applied to the imaging 





of stationary water-borne targets. The results of computer 








simulations designed to compare different methods of esti- 


mating acoustic signal propagation ranges and received 








phase angles are presented in Chapter IV. Finally, Chapter 
V presents the thesis conclusions and recommendations for 


future research in SAS signal processing. 








The next chapter describes the treatment of the ocean 





as a linear, time-variant, space-variant filter. The SAS 
mine-hunting geometry and scattering model are also intro- 


duced. 





TH 


IS PAGE 








NTENT 





ONALLY LEFT BLANK 





II. THE OCEAN AS A LINEAR, TIME-VARIANT, SPACE- 
VARIANT FILTER 


This chapter presents an overview of the application 


of linear systems theory to acoustic propagation modeling. 


A. INTRODUCTION 
The computer simulations presented in this thesis are 
based on a linear systems model of the ocean environment. 


A linear systems approach to acoustic pulse propagation is 





presented in [4] and [5]. All simulation models in this 








thesis are based on treatment of the ocean as a linear 


filter. 


Propagation of small amplitude acoustic pulses in the 





ocean can be described by a linear wave equation. This 





linear wave equation accurately describes the propagation 


of acoustic pulses from the source to the discrete point 





scatterer and from the discrete point scatterer to the 


receiver. To model this propagation, we treat the ocean as 





a linear, time-variant, space-variant filter. 


B. COMPLEX FREQUENCY RESPONSE OF THE OCEAN 


The propagation of small amplitude acoustic pulses 





through the ocean medium can be described by the linear, 





three-dimensional, inhomogeneous wave equation, [5] 





oO” 


1 
Vi lEN)— Sey ap Iu EN = tulen) s (D515 


where x,(t,r) is the source distribution at time ¢ and posi- 
tion r with units of inverse seconds and jy,(¢t,r)is the veloc- 


ity potential of the acoustic field at time ¢ and position 





r with units of squared meters per second. The vector 
r=(x,y,z) describes the location of any point in three- 


dimensional space. The position dependent speed of sound 








in meters per second is described by c(r). Note that both 
the source distribution and the velocity potential are 
functions of time and position. The time-variant property 
of the linear filter allows us to model motion of the 
transmitter, target, and receiver. This ability to accu- 


rately model motion of all platforms is essential in SAS 








modeling where transmitter and receiver motion is used to 








synthesize the array. This time-variant property also is 


used to account for changes in the ocean medium with time. 





The space-variant property allows for the presence of 





boundaries, discrete point scatterers and spatial variation 


in ambient density and the speed of sound. 


The solution of Eq. (2.1) can be obtained by treating 





the ocean as a linear, time-variant, space-variant filter. 


Figure 4 is the linear system block diagram of the ocean as 





a time-variant, space-variant filter. The input-output 


relationship is given by [5] 





yy (tr) = ec ar (tr |t,.1, )dt,dr, , (2,2) 


where hg lielin jas the time-variant, space-variant, impulse 
response (Green’s function) of the system. The impulse 


response fixe (it, |t5.%q.) ie “the Tesponse OF “Lhe: Tiler at eime.t 
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and position r=(x,y,z) due to application of a unit ampli- 
tude impulse at time 4% and position ry =(%).¥%,%)-. The input 
acoustic signal x,(t,r) is the source distribution at time t¢ 


and position r and jy,(¢,r) is output acoustic signal at time 





t and position r. It is important to realize that Eq. (2.2) 


is a four-fold integral since dr, =dx, dy, dz,. 


hy (t.r|t).%,) 


xy(t0) Vy (40) 


H,,(t.r| f.v) 





Figure 4. A linear, time-variant, space-variant, ocean 
medium filter [From Ref. 5]. 


Figure 4 also shows the time-variant, space-variant, 
transfer function of the ocean filter H,,(t.r|f,v) defined by 


[5] 
foie (tr |f.v) = | a (tr | f.% )exp[+j2zv © (r—1,)] dr, 7 (2.3) 


where f represents input frequency components in Hertz and 


v=(f,.f).f,) is a vector representation of the spatial fre- 





quencies in cycles per meter in the xX, Y, and Z directions, 
respectively. The time-variant, space-variant, complex 


frequency response of the ocean at frequency f is [5] 


Hy (tr | fot) = [Au (tr | foot) expl—j2af(t-t)l at, « (2.4) 
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The strategy fol 





the time harmonic sol 


lowed in Reference 5 is to first find 


ution to the linear, inhomogeneous, 





three-dimensional wav 





tion is then used to 

equation by using Fou 
The focus of thi 

response of the ocean 


Eq. (2.4). The comple 


medium integrates th 





path between transmit 





e equation. The time harmonic solu- 





find a pulse solution to the wave 


rier transform techniques. 





s thesis is on the complex frequency 
medium, Biel ta) as defined by 
x frequency response of the ocean 


characteristics of the propagation 


ter, discrete point scatterer, and 





receiver. Simulation 
approximate propagati 
gation delay as deter 


acoustic field. 


CG. SUMMARY 

The ocean medium 
variant, space-varian 
develop a time-varian 


response of the ocean 





signal propagation an 


linear, inhomogeneous 





C22) In the next c 


s were conducted using both exact and 


on ranges to account for signal propa- 





mined by propagation path of the 


can be modeled as a linear, time- 


t, filter. This approach allows us to 





t, space-variant, complex frequency 


that accurately describes acoustic 





d therefore provides a solution to the 





, wave equation given by equation 


hapter, the complex frequency response 


of the ocean is derived for the specific case of SAS ap- 


plied to stationary o 


bject detection. 
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III.MODELING OF BISTATIC SCATTERING FOR SYNTHETIC 
APERTURE SONAR 


This chapter introduces the mathematical model used 
for simulation of synthetic aperture sonar (SAS) as applied 


to the imaging of stationary water-borne targets. Several 








traditional approximations and an exact solution for propa- 
gation ranges are introduced and the limitations of each 


are discussed. 


A. INTRODUCTION 


The main advantage of a SAS system is the ability to 





synthesize an extremely long array while maintaining a 





relatively small actual size of the physical array. Syn- 





thetic aperture sonar beamforming and signal processing 





relies on properly steering and focusing the aperture beam 





pattern in order to co-phase all the received signals. Due 








to the effects of motion in the SAS problem, the propaga- 


tion path from the transmitter to the discrete point scat- 











terer and from discrete point scatterer to the receiver is 


time varying. Traditionally, simple approximations are 











used to determine these propagation ranges and angles of 





incidence and scatter. This chapter introduces the exact 





bistatic scattering model for SAS systems presented in [9] 
as well as some common approximations used when modeling 


acoustic signal propagation in SAS systems. 
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B. BISTATIC SCATTERING MODEL 


1. Assumptions 


The focus of this section is to develop 





frequency response of the ocean for a SAS system with 











speed of sound and the ambient density of 


treated as constants. This allows 


tionary targets. In the binomial 








the exact propagation models described in this thesis 





the complex 


sta- 


approximation models and 


, the 


the ocean are 


us to treat the ocean 


medium as homogeneous. AS a result, sound rays will travel 








in straight lines from the transmitter to the discrete 

















point scatterer and from the discrete point scatterer 





the receiver. Propagation model 





EO 


ls involving surface and 


bottom interactions can be developed in an identical manner 





to those presented for propagation between the 





transmitter 


and discrete point scatterer and discrete point scatterer 


and receiver [9]. 


The transmitter and receiver velocity vectors are as-— 








sumed to be constant only during the transmission and re- 


ception of the acoustic pulse. 


The platform may change its 


speed at any time in any direction between transmission and 


reception and between cycles. 


need for motion compensat 


signal processing algorit 





thms. 


This feature eliminates the 


tion typically applied to most 





The binomial approximat 


t SAS 





tion 





and exact propagation models assume only that the position 


and velocity vectors of 





determined during signal transmission and reception. 





termination of the vehic]l Vv 


the transmitter and receiver can be 





De- 


locity vector is not a trivial 





problem, however, and accurate det 
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termination of vehicle 





position is an area of extensive research in SAS signal 


processing and vehicle design [1, 10, 11]. 








In order to allow comparison of the various methods of 
propagation range determination, the scattering function of 


. . A AN . 
the discrete point scatterer, elf mene) is assumed to be 


unity. 


2. Propagation Geometry 
We begin our development of the complex frequency re- 


sponse of the ocean by defining the velocity vectors for 


the transmitter, V,, the diserétée point searterer, V,, and 








the receiver, V,. These velocity vectors are given by [9] 
Vy =Vomy, + (3.1) 
V,=Vin, + (3.2) 
and 
Vi =Viny, + (693) 


where Vo>. V9 and VY. are the speeds in. meters per second of 





the transmitter, discrete point scatterer, and receiver, 





respectively. The dimensionless unit vectors fy,, my, and 





mn, define the directions of V,, V,, and V,, respectively. 


Velocity vectors given by Egs. (3.1) through (3.3) are con- 














stant as discussed in Section -B.1. Motion is consid- 





ered to begin at time f=f, . 


The time-varying geometry of the scattering model is 





modified using relative velocity vectors to allow the 


transmitter and the discrete point scatterer to be treated 
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as being motionless. The velocity vector of the discrete 





point scatterer relative to the velocity vector of the 
transmitter in the direction of the velocity vector of the 


discrete point scatterer Ny is [9] 


Vio =Vi— (i, © Vo) ty, =| Vo (iy 9 ty, ) |i, - (3.4) 


i 


The initial propagation model allows for motion of the 





discrete point scatterer, so we define a relative velocity 





vector, Vis The vector Vy, 1s: the velocity vector of the 


receiver relative to the velocity vector of the discrete 
point scatterer in the direction of the velocity vector of 


the receiver ny, and is given by [9] 
V2, = V2— (ty, °Vi) a, =[V2-V, (4, 9%) Ay, - (3.5) 


The position vectors from the origin to the transmit-— 


ter, discrete point scatterer, and receiver when motion 





begins at t=t, are given by rm =(%).M.2%), N=(%>2%), and 


r, =(X,,,,Z,), respectively. Using these position vectors, we 





can further define the position vector from the transmitter 


to the discrete point scatterer as [9] 
mH =n-h (3.6) 


and the position vector from the discrete point scatterer 


to the receiver as [9] 
r.=r-7,.- boat) 


The dimensionless unit vectors from source to discrete 





point scatterer and from discrete point scatterer to re- 


ceiver are [9] 
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hie (356) 
ro,| 
and 
: r 
a (3.9) 
In,» 
respectively. The ranges in meters from the transmitter to 


discrete point scatterer and from the discrete point scat-— 








terer to receiver are given by [9] 





7 = Nou (3.10) 
and 

ho Ne ’ (3.11) 
respectively. This bistatic scattering geometry is shown 


in Figure 5. 








It is now important to define the parameters of our 





model with respect to the time instance when the acoustic 





field is first incident upon the discrete point scatterer 





and the time instance when the acoustic field is first 





incident upon the receiver. The transmitted acoustic field 








is first incident upon the discrete point scatterer at some 





time ¢t/ seconds where ¢’/>t,. This allows us to define the 





time difference from transmission of the acoustic field to 





first incidence upon the discrete point scatterer as [9] 
Meta. £ >t (2212) 


Similarly, the time instant when the scattered acoustic 


field is first incident upon the receiver is t>f>t,. The 
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time difference between ¢ and first incidence of the acous- 





tic field upon the receiver is given by [9] 


At’=t-t, t>f>t,. (3.13) 





Figure 5. Bistatic scattering geometry when motion begins 
at time t=t, seconds. Point 0, P(r), is the transmitter; 


point 1, P(r), is the discrete point scatterer; and point 


2, P(r), is the receiver. All three platforms are in mo- 
con. (From Ret, 5]. 


By referring to Figure 6 we can express the position 


vector from the point source to the discrete point scat-— 


terer at time f as rj, and the position vector from the 
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discrete point scatterer to the receiver at time ¢ as r,. 





The position vector from the origin to the discrete point 





scatterer at time f is r. The: positon: vector from the 


‘ ‘ ‘ : F; ” 
origin to the receiver at time f is r, . 


3. Velocity Potential and Complex Frequency Response 
Our goal is to develop an expression for the complex 


frequency response of the ocean in order to compare differ- 





ent methods for estimation of propagation range from trans-— 








mitter to discrete point scatterer and from discrete point 





scatterer to receiver. Prediction of the propagation 
ranges Ir, and Ir’, | will allow us to correctly steer and 


focus the beam pattern of our synthetic array resulting in 











the co-phasing of the output electrical signals at the 











different receiver locations. Once a model for the complex 





frequency response is developed, we can predict the acous- 





tic field (velocity potential) at the receiver. 


When the transmitted acoustic field is first incident 











upon the discrete point scatterer at time f’, the position 





vector from the transmitter to the discrete point scatterer 


is given by [9] 


ro. =Mo t At’ Vip - (3.24) 
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S,0(r - 4) 


S,0(r - 5) 





Figure 6. Bistatic scattering geometry when the transmit-— 
ted acoustic field is first incident upon the discrete 
point scatterer at time ft seconds and when the scattered 
acoustic field is first incident upon the receiver at time ft 








seconds where ¢t>f>t,. Point 0, P(r), is the transmitter; 
point 1, P(r), is the discrete point scatterer, and point 


2, P(r), is the receiver. All three platforms are in mo- 
tion [Prom Ref. 5] « 
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Recall that we are modeling the propagation of the acoustic 
field through the ocean as transmission through a linear 


time-variant, space-invariant filter. Therefore, th 





acoustic field incident upon the discrete point scatterer 





at time ¢’ and position Slee) is [9] 


Vu (*.1/) Sanh sie Ea |/ )exp(+j27 7’), be (3.15) 
where 


pa ieleal Bites 


7 vd 
Felina : 
4z|r,,| 





is the time-variant, space-invariant, complex frequency 


response of the ocean at frequency f hertz. The source 


strength S, is given in cubic meters per second and k is 











the wavenumber in radians per meter given by 
R= 27 f (e=27) Ax (S27) 


An identical method is used in [9] to develop an ex- 


pression for the acoustic field incident upon the receiver 











at time ft. The position vector between the discrete point 





scatterer and the receiver at time oe aes is given by [9] 
r, =r, + At’ (V,-V,.)+AL’V,, . (3.18) 


The acoustic field (velocity potential) incident upon the 


receiver at time ¢>f>t, is [9] 


Vij (1.1,")=S\Hty (1,13 |f.r,)exp(+j2z ft), t>t >t, (3.19) 


where 
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de ioe (1,17 fh aH, (4, Sjedineno) ty (6% If) 
ay Ox | — Fe (al +fe%) eee) 
= & (f.7%1.45) 162 r, ry, . 














The complex scattering function of the discrete point 
scatterer, g,(f,%%,,4,), is dependent upon frequency and upon 


the angles of incidence and scatter. For simplicity, this 





function is assumed to be unity in this thesis so that we 





may independently investigate the phase term in Eq. (3.20). 


Cc. PROPAGATION RANGE DETERMINATION 


The key to proper co-phasing of the output electrical 





signals from the SAS array is estimation of the received 


phase of the acoustic field incident upon the receiver. 





Examination of Eq. (3.20) reveals that the phase term, 





exp| — jk (|ra,|-+[r7a]) |, depends directly on frequency and the 





magnitude of the position vectors rj, and rj,. It is in the 





solution for these time-varying ranges that this thesis is 








focused. Accurate determination of the propagation range 


of the acoustic signal is especially important in recent 











SAS projects such as the Long-Term Mine Reconnaissance 








System (LMRS) where relatively high carrier frequencies 
will be used to increase image resolution. Small errors in 


propagation range estimation can introduce large errors in 








signal phase estimation at the receiver, significantly 


degrading the array performance. 
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1. The Stop-and-Hop Approximation and the Moving Re- 
ceiver Correction 


One approach to estimation of propagation ranges is to 








assume that the transmitter and receiver are stationary 


during signal transmission and reception. This assumption 





is acceptable for many synthetic aperture radar (SAR) sys- 





tems, but due to the slow propagation speed of acoustic 


waves in water relative to electromagnetic waves in air, 








this assumption is invalid for a SAS system imaging targets 


at any significant range [6]. An initial correction may be 





applied by compensating for the motion of the receiver in 











the time interval between signal transmission and signal 
reception. A method for received phase correction due to a 


moving receiver is presented in [6]. Figure 7 shows the 





geometry involved in this correction. The initial range 





from transmitter to targét is.-K,. The initiel range from 





target to recéiver is: Ky. The distance traveled by the 


receiver during signal propagation is Vt) where V is the 


vehicle speed and ¢’ is the total propagation time between 





transmission and reception of the acoustic signal. The 


propagation time ¢ is an estimation of the total propaga- 








tion time from transmitter to target and from target to 





receiver based on a reference range to the scene center. 


The range R, describes the range from target to receiver at 


time ¢/. Each received acoustic signal is multiplied by the 





phase correction term [6] 


exp( 4/2), (5.20%) 
Cc 
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Be 
Receiver 


vy 


Figure 7. The moving receiver correction for the stop- 
and-hop geometry [After Ref. 6]. 


where 


AR=R,-R, . (3.22) 


Compensating for the motion of the receiver during 





signal propagation does remove some error inherent in the 


stop-and-hop assumption but several additional factors are 





not considered. The motion of the transmitter and receiver 





during signal transmission and reception is ignored. In 








addition, the vehicle velocity vector is assumed to be 
constant for the entire duration of signal propagation. 
This correction further assumes that the total propagation 
time, ¢ in Figure 7, from transmitter to target and target 


to receiver can be approximated by calculating the propaga- 





tion time to some reference range, R,, such that [6] 


t=—_. (3.23) 





The reference range, K,,. is the range to the center of 


the image area of the SAS system [6]. Throughout the re- 








mainder of this thesis, the moving receiver correction is 














applied to the stop-and-hop model and referred to as the 


stop-and-hop approximation model. 


We can define the stop-and-hop model presented in [6] 

















in terms of the bistatic scattering model presented in 

















Section -B of this thesis. The range from transmitter 








to discrete point scatterer is given by 


R, = Ir, -h 





, (3.24) 


Where © <and 1, arevthe positron sectors Erom the Goordinatre 








system origin to the discrete point scatterer and transmit-— 


ter, respectively (see Fig. 5). The time-varying range 








from discrete point scatterer to receiver is given by 





ae, a ae (3.25) 





R,(t)=[r, (1-4, 
where r,(t) is the time-varying range from the coordinate 
system origin to the receiver given by 

nia=an()t+Vt, >. (3.26) 


The vector V is the velocity vector of the receiver, and ¢’ 


is given by Eq. (3.23). 


2. The Binomial Approximation Method 


The magnitude of the position vectors rj, and rj, can be 





approximated by using the binomial expansion presented in 


[5]. The magnitude of the position vector from the trans-— 
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mitter to the discrete point scatterer at time ¢f can be 


rewritten as [5] 








pe =nVivb =n, (145-24-] jo]<1 (3597) 
where 
VioAb ise aoe A 1 VAL 
b=2—"—_| (fi, 7, )+-—"— | . (3.28) 
im : pet M1 


Recall that in our bistatic scattering model, t is the time 








instant when the transmitted acoustic field is first inci- 





dent upon the discrete point scatterer and that Ute For 





the case of a SAS system imaging a stationary target, 
Vi,=0. (3.29) 
Therefore, Via given by Bq. (3.4) 2s undefined, As a result, 
Vio =9 (3.30) 
and 


Viol =o 20's (3.31) 





In this case, b given by Eq. (3.28) is equal to zero and 


Eq. (3.27) reduces to 


, 


Ny 








=r) (3.32) 


a 


which is exact. 





The binomial approximation can also be applied to the 


estimation of Ir’, | « The range Ir’, | can be written as [5] 
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where 


2 


b=] fi, °(V.-V,o) ]At’+(7,.° V2) 7} . (3. 


No 


The expression for b can be simplified using 
Eqs. (3.29) and (3.30) for the stationary target. Thus, 


Eq. (3.5) reduces to 


V,,=V; (3. 


and Eq. (3.34) can be rewritten as 


b=—{(Ar'+Ar’)(i,, ev,)}. (3. 


No 


Using Eqs.(3.12) and (3.13), 


At +At’=t-t,+t-f =t-t, =At. (2% 


The expression for b now reduces to 


2 


b=—{Ar(A,,°V,)}. (3. 


No 


Equation’ (3.33)) as-only Valid Lt jb|<1. Note that 





la. °V.|<|Vo] - (3. 


Now, if we impose the more stringent criterion that 


2/V,|At — 


Ko 


a Be (3. 


é3:3:) 


34) 


3:5.) 


36) 


37) 


38) 


39) 


40) 


then the binomial approximation can be simplified by using 











only the first two terms of Eq. (3.33). Doing so yields 
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ra |= 7.2 +] An °(V2- Vio) Ae +(a,.°V,,) Ae” (3.41) 


Equation (3.41) can be further simplified using Eqs. (3.30), 
(3.35), and (3.37) such that 


=1,+(fi,)°V,) At . (342) 








32 Exact Propagation Range Method 


Reference 9 presents a method for determining the 


time-varying propagation ranges Ire (0) and r7,(0| exactly for 





the general case of bistatic scattering when all three 





platforms (transmitter, target, and receiver) are in motion 
and for the special case of a SAS system imaging a station- 


ary target. The constant value of range from transmitter 











to discrete point scatterer when the transmitted acoustic 








field is first incident upon the discrete point scatterer 





at time instant >t, is given by [9] 


, 


Yo) 








=1,. (3.43) 


The constant value of range between discrete point scat-— 


terer and the receiver when the scattered field is first 











incident upon the receiver at time instant ft>?f’>t, is given 








by [9] 
yi Bot B, F4aC, 
(SS (3.44) 
; 2A, 
where 
2 
eas, (3.45) 
(6 


28 





























7 eV y2 
B22 oe ae Don, oy (3.46) 
“C 
and 
Vo (7, © Vo) 
Co= Mo crt Mouths — ue Vea24 ty 
The time-varying propagation ranges are given by [9] 
” 2 a 2 27 1/2 
IrOl=[".+2n2G,eVAttM (Ad |, t2t,4+7, (3.48) 
and 
[ri] = cAt-]r70) » €24,+7, (3.49) 
where 
r r, 
me 0,1 re 1,2 (3.50) 
Cc Cc 


is the time delay in seconds or the 


for the transmitted acoustic signal 








the receiver after motion begins at 


At=t-t,, t2t,+7. 


D. SUMMARY 


time instant ¢,, 


amount of time it takes 


to begin to appear at 


and 


(3.51) 


Acc 


urate estima 


tion of the propagat 


tion range from the 








transmit 


ter to targe 


t and from target to receiver is essen- 














tial to 





signals from the different receiver locations. 








tional stop-and-hop method, 


can be used to estimate propagation ranges. Al 





the propagation ranges can be determined using 
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the proper co-phasing of the output electrical 


The tradi- 


corrected for receiver motion, 


ternately, 





the bistatic 


scattering model presented in this chapter. The propaga- 





r/,(0| can be estimated using a bino- 





tion ranges Ir, ()| and 


mial approximation or calculated exactly. The next chapter 











presents the results of computer simulations used to com- 
pare the accuracy of the stop-and-hop approximation, bino- 


mial approximation, and the exact propagation range calcu- 


lation methods. 
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IV. PROPAGATION RANGE AND RECEIVED PHASE ANGLE 
ESTIMATION COMPARISON 


This chapter presents the results of computer simula- 
tions designed to compare different methods of estimating 


acoustic Signal propagation ranges and received phase an- 





gles for a SAS system imaging a stationary target. Simula- 


tion results of the receiver output electrical signals for 





a SAS array with Fast Fourier Transform (FFT) beamforming 


and various propagation range estimation methods are also 


presented. 
A. INTRODUCTION 
1: Method of Evaluation 


Accurate estimation of the propagation range from the 


transmitter to the discrete point scatterer and from the 





discrete point scatterer to the receiver is essential to 








properly co-phase the output electrical signals in the SAS 














system. Using the models presented in Chapter , the 





simulations presented in this chapter compare the propaga- 
tion ranges and received phase angles computed using the 

binomial approximation method and stop-and-hop approxima- 
tion method to the exact propagation ranges and exact re- 


ceived phase angles. Comparisons are based on the percent 








error in range from transmitter to discrete point scat-— 

















terer, percent error in range from discrete point scatterer 





to receiver, and percent error in the estimation of the 





phase of the signal incident upon the receiver. All simu- 


lations were conducted using MATLAB version 6.0. 
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Comparisons of propagation range and received phase 


estimation methods are performed for a single trans-— 








mit/receive cycle. The propagation range error for the 











time-varying range from transmitter to discrete point scat-— 





terer is given by 





ri) 








ri). (4.1) 


Exact Approx 


Similarly, the propagation range error for the time-varying 


range from discrete point scatterer to receiver is given by 


IA) 





A ae (4.2) 


Exact Approx 


The percent range error for the time-varying range from 





transmitter to discrete point scatterer is given by 








Mou Var s Mou (0) 


: “Pre* %100% . (4.3) 
mi O 











Exact 


The percent range error for the time-varying range between 


discrete point scatterer and receiver is given by 





1, (t) 








” 
Exact - Me (| 


I, (¢) 





APPEL IQO% . (4.4) 





Exact 





Finally, the phase in radians of the acoustic signal inci- 


dent upon the receiver is given by 


a(1)= =| 








(0+ |r7.(0)) (4.5) 


and the percent phase error is given by 


P(t) rracr ~ (1) 


ApPre* 100% . (4.6) 
0 (isin 
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2. Simulation Parameters 
The computer simulations developed for this thesis 


simulate vehicle motion in any direction in an unbounded, 





homogeneous, ocean medium. Although the simulation places 

















no restrictions on target location, the target is always 








assumed to be at some depth below the transmit/receive 


platform. Unless otherwise noted, the vehicle begins mo- 








tion. ati time. 2. ab: fhe origin ofthe -coordinate system: and 


is imaging a target forward of the beam (@=45°) at a range, 


R (see Fig. 8). For simplicity, platform motion is along 





the positive Z axis. 


(r.@.¥) 





Figure 8. SAS simulation geometry [After Ref. 5]. 
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Propagation ranges are calculated with the equations 

















presented in .C.1 through .C.3 using the assumptions 























presented in +B 





The platform speeds and target ranges used for com- 





parison are similar to those that have been specified for 
the Navy’s LMRS system [2]. The LMRS side-scan sonar sys- 
tem is used to synthesize the SAS array. The side-scan 
sonar array covers both sides of the six-foot long LMRS 


platform, allowing independent imaging on both the port and 


starboard sides of the vehicle. The vehicle is designed to 
travel between 3 and 7 knots [2]. The minimum speed is 
based on maintaining maneuverability. The maximum speed is 





set to ensure adequate sampling of the target area. Figure 
9 shows the relationship between physical array length, 


platform speed, and maximum range. 


MAXIMUM SAS COVERAGE (PER SIDE) FOR VARIOUS 
PLATFORMSPEEDS AND ARRAY LENGTHS 


HIGH-SPEED PLATFORMS RE QUIRE 
20 LONGER PHYSICAL ARRAYS TO 
ACHIE VE LONG-RANGE COVERAGE 


_ 
“ 


10 


LMRS ARRAY LENGTH 


SPEED (KNOTS) 





0 200 400 600 800 1000 
RANGE (YARDS) 
Figure 9. LMRS Range and Speed limitations [From Ref. 2]. 
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3. FFT Beamforming 

The simulations presented in the next section compare 
the accuracy of the approximate stop-and-hop propagation 
range method and the approximate binomial approximation 


propagation range method for each transmit/receive cycle of 








the platform. To provide further comparison and to better 





analyze the effects of approximation errors on SAS systems, 





the propagation range calculations from the stop-and-hop 


approximation, binomial approximation, and exact methods 














are used to estimate the phase weights in a Fast Fourier 





Transform (FFT) beamforming algorithm. 


To independently investigate the performance of each 





method due to phase errors, we assume the scattering func-— 





tion of the target is unity as discussed in Sections 
































.B.1 and .B.3. We also neglect the propagation range 


term, l6a’|rj,||r7,{, in the denominator of Eq. (3.20). Neglect- 














ing this amplitude factor allows us to focus on the effect 


of propagation range error on received phase prediction. 





Three simple sinusoids, each at a different frequency, 





are used to simulate acoustic propagation. As a result, 


Eq. (3.19) can be expressed as 








3 
vu (05 (t)) = Vexp| = sk; (|r (0|+[r2.(0)) | exp(i22- fe), (4.7) 
i=] 
where 
k, Bll (4.8) 
Cc 


is the wavenumber for each different frequency and f,=10 





kHz, f,=20 kHz, and f,=30 kHz. 
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We assume that the output electrical signal of the re- 








ceiver is directly proportional to the acoustic field inci- 








dent upon the receiver. Therefore, 


(452) =u (tr7() . (4.9) 


We now describe the FFT beamforming algorithm used in 





this thesis and is shown in Fig. 10. FFT beamforming is 
accomplished by first taking the time-domain Fourier trans-— 


form of the output electrical signals at the receiver at 

















different positions as the SAS system images a target at a 








known location. These signals are simulated using the 


exact propagation model to compute the acoustic signal 


incident upon the array given by Eq. (4.7), where Ir, ()| and 


In, are given by Eqs. (3.49) and (3.48), respectively. The 





resulting output electrical signals are given by Eq. (4.9). 


These exact electrical signals are then phase weighted with 





phase weights computed using one of the approximate meth- 


ods. 


These approximate phase weights are computed by first 





determining the acoustic field incident upon the receiver 





using Eq. (4.7), where the values for Ire (0)| and AG) are 





computed using either the binomial approximation method or 





the stop-and hop approximation method. The approximate 
stop-and-hop values for Ir, (| and A) are given by 
Eqs. (3.24) and (3.25), respectively. The approximate bino- 


mial approximation values for Ire (0)| and AG) are given by 





Eqs. (3.32) and (3.42), respectively. The negative of the 


phase spectrum of the time-domain Fourier transform of the 
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output electrical signals given by Eq. (4.9) determine th 





appropriate approximate phase weights to use to phase 


weight the receiver output electrical signals. 








T/R Cycle 1 


T/R Cycle 2 














T/R Cycle N 








Figure 10. FFT Beamforming Algorithm. 


The approximate phase weights are then multiplied by 


the time-domain Fourier transform of the output electrical 
Oo 


signals determined using the exact method. The exact out- 





put electrical signals have now been phase weighted using 





approximate phase weights. We then compute the invers 





Fourier transform of the phase-weighted signals. The re- 
sulting time-domain signals are added to produce the total 


output electrical signal from the SAS array given by 





N 
yr(t)= Dye, (t-(t,, +%)) + (4.10) 
n=l 


where N is the total number of transmit/receive cycles used 


to create the SAS array and Ye, (t-(t, +,)) 16 Ene Out pue 


electrical signal from the receiver for the n” transmit/ 








receive cycle after phase weighting. The signals are syn- 





chronized in time by subtracting the time delay for each 


transmit/receive cycle, ti +,), where ¢, is the time at the 





beginning of the n” eyelé-and: 7 Ts given Dy Bos(3.50): In 


the FFT beamforming algorithm shown in Fig. 10, we assume 








that the output electrical signals have been synchronized 











in time before performing the forward FFT. With the trans-— 


mit/receive platform constantly in motion, the initial 





range from the transmitter to the discrete point scatterer 





and from the scatterer to the receiver will be different at 





the beginning of each transmit/receive cycle. These ini- 








tial ranges are updated to reflect the geometry at the 
beginning of each cycle and, as a result, the propagation 


delay, T,, is recomputed for each transmit/receive cycle. 


The FFT beamforming procedure is described in detail in 


[4]. 
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B. PROPAGATION RANGE AND RECEIVED PHASE COMPARISON 


1. Stop-and-Hop With Moving Receiver Correction 


The performance of the stop-and-hop propagation model 

















is evaluated using Eqs. (4.1) through (4.6) where Ire], _ is 

given by Eq. (3.49), ri (0)|, is given by Eq. (3.24), [LO 
> [pprox ? xaci 

is given by Eq. (3.48), and AG) is given by Eq. (3.25). 


Approx 


Figure 11 shows the errors introduced by the stop-and- 


hop approximation for a vehicle speed, V,, of 2-7 knots, 





R=213 yards, and frequency, f, of 100 kHz for a single 





transmit/receive cycle. In this case, the stop-and-hop 








approximation with the moving receiver correction intro- 


duces no error in the estimation of the range from discrete 





point scatterer to receiver, In”, O] - The error in the esti- 


mation of range from transmitter to discrete point scat- 





terer, Ta) 





, is, however, significant. The error intro- 


duced by the stop-and-hop method results in a maximum per- 





cent received phase error of -0.5% at V,=7 knots. 


At the minimum platform speed of 3 knots, we can in- 





vestigate estimation errors from the stop-and-hop approxi- 





mation over all possible target ranges observing the limi- 








tations described in Figure 9. Figure 12 shows propagation 
range and received phase error for V,=3 knots, 


f =100 kHz, with R varied from 25-450 yards. From Figure 





12 we can see that the estimation error for Ir, (| and the 





subsequent phase error increases as the range to the target 
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decreases. This suggests that the stop-and-hop approxima-— 
tion provides unacceptable performance inside some minimum 


range. 


Ir 4)I Error: Stop—and—Hop Io MI % Error: Stop—and—Hop 
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Figure 11. Propagation range and received phase error 
using the stop-and-hop approximation for V,=2-7 knots, 


f =100 kHz, and R=213 yards. 
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Figure 12. Propagation range and received phase error 
using the stop-and-hop approximation for V,=3 knots, 
f =100 kHz, and R=25-—450 yards. 
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It is interesting to examine the accuracy of the stop- 





and-hop approximation with the target at different initial 








angles with respect to the array. The angle @ is a cylin- 


drical coordinate measured from array broadside to the 








target projection in the X-Z plane at the time motion be- 





gins, as seen in Figure 8. This simulation is conducted 


such that R remains constant at 300 yards while @ is var- 





ied from 0° to 90°. The target depth is set at 150 yards 


below the platform. ithe platiorm Speeds. Voy 18) S -KNOES « 





Figure 13 shows that propagation range and received phase 
error increase as @ increases. This is an important obser- 
vation since all SAS systems employ the coherent addition 
of multi-aspect data to form high-resolution images. The 


ability to accurately image a target at higher aspect an- 








gles results in improved along-track image resolution [2]. 
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Figure 13. Propagation range and received phase error 
using the stop-and-hop approximation. Angle @ varies from 


0° to 90°. Speed V,=5 knots, f=100 kHz, and R=300 yards. 
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2. Binomial Approximation Method 


The performance of the binomial approximation method 








for propagation range estimation is evaluated using Eqs. 


(4.1) through (4.6) where raO),.. is given by Eq. (3.49), 





Ba Ol.. is given by Eq. 





ri()|, is given by Eq. (3.32), 
> [pprox 
(3.48), and 2 OL pos is given by Eq. (3.42). 


Simulations examining the accuracy of the binomial ap- 


proximation model were conducted in the same manner as 





those for the approximate stop-and-hop model. Once again, 








the binomial approximation estimates of the propagation 





ranges |r (t)) and |r7,(f| and the resulting estimated phase of 
0,1 12 





the acoustic signal incident upon the receiver were com 





pared with the exact solutions for these parameters. 


Figure 14 shows the errors introduced by the binomial 





approximation for a vehicle speed, V,, of 2-7 knots, 


R=213 yards, and frequency, f, of 100 kHz for a single 





transmit/receive cycle. In this case, the binomial ap- 








proximation introduces extremely small error in the estima-— 


tion of the range from the transmitter to the discrete 





point scatterer, Ir, (0) . Maximum estimation errors for the 





binomial approximation of Ir, (| are on the order of 10° 








meters at =) knots. The error from. the estimation oF 
range from discrete point scatterer to receiver, Ir.) is 
also extremely small. Maximum estimation errors for the 


binomial approximation of Ir’. are on the order of 10° 
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meters at V,=7 knots. The resulting phase error introduced 


by the binomial approximation method results in a maximum 





percent received phase error of 0.004% at V,=7 knots. 





Recalling that the stop-and-hop approximation resulted ina 


maximum percent phase error of —-0.5% for the same scenario, 








it is obvious that the binomial approximation greatly re- 
duces propagation range and received phase estimation er- 


EOE. 


At the minimum platform speed of 3 knots, we can again 





investigate estimation errors from the binomial approxima-— 


tion over all possible target ranges observing the limita- 





tions described in Figure 9. Figure 15 shows propagation 


range and received phase error for V,=3 knots and f=100 








kHz while R is varied from 25-450 yards. From Figure 15 we 





can see that the estimation error for Ir, (| increases as the 


target range, R, increases. Conversely, the range estima- 





ELON “Sr ror -for AG) decreases as R increases, and is much 


greater than the range estimation error for Ir, (| Lp Pala 


cases. This results in a phase error that increases as R 











decreases, a result similar to that observed for the stop- 
and-hop approximation. This suggests that the binomial 


approximation also may provide unacceptable performance 








inside some minimum range. Comparing the stop-and-hop 


approximation method with the binomial approximation 














method, we again see that the binomial approximation method 








introduces much less error in the estimation of propagation 


ranges and received phase. 
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Figure 14. Propagation range and received phase error 
using the binomial approximation method for V,=2-7 knots, 


f =100 kHz, and R=213 yards. 
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Figure 15. Propagation range and received phase error 
using the binomial approximation method for V¥,=3 knots, 


f =100 kHz, and R=25-—450 yards. 
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We can also examine the accuracy of the binomial ap- 





proximation method with the target at different initial 





angles with respect to the array. The angle @ is a cylin- 


drical coordinate measured from array broadside to the 








target projection in the X-Z plane, as seen in Figure 8. 


This simulation is conducted such that R remains constant 








at 300 yards while @ is varied from 0° to 90°. The target 





depth is set at 150 yards below the platform. The platform 
speed, V,=5 knots. Figure 16 shows that propagation range 
and received phase error decrease as @ increases. This 


suggests that the binomial approximation method would be 


able to more accurately combine the multi-aspect data in- 





herent to SAS systems. Once again, we can see that the 





binomial approximation introduces a maximum of 0.002% phase 





error compared to -0.25% phase estimation error per trans-— 





mit/receive cycle introduced by the stop-and-hop approxima- 


tion. 
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Figure 16. Propagation range and received phase error 
using the binomial approximation. Angle @ varies from 0° 


to 90°. Speed V,=5 knots, f=100 kHz, and R=300 yards. 
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Ci FFT BEAMFORMING SIMULATIONS 
To provide some insight into the performance of the 
various propagation models when applied to a SAS system, 


each model was used to calculate the phase weights for use 





in an FFT beamforming algorithm as described in IV.A.3, 


where f,=10 kHz, f,=20 kHz, and f,=30 kHz. 


The array is synthesized as the vehicle travels at the 


designated speed in a straight line along the positive Z 





axis. 


Figure 17 shows results of FFT beamforming with 


(a) Vj=3 knots and (b) Vj=7 knots. The array is synthe- 








sized using 10 transmit/receive cycles with the target at 
an initial range, R, of 213 yards and @=45°. Phase weight- 


ing the output electrical signals using the binomial ap- 














proximation method to determine received phase of the inci- 
dent acoustic signal at both 3 and 7 knots produces a wave- 
form that is almost exact. The errors in the stop-and-hop 


approximation, however, are significant and phase weighting 





the output electrical signals using this method results in 





very low array gain. The poor performance of the stop-and- 


hop approximation results from the relatively high error in 





the received acoustic signal phase estimate. 
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Figure 17. FFT beamforming with (a) V,=3 knots and (b) 


V,=7 knots. 


receive cycles with the target at an initial range, 
213 yards and @=45’. 


The array is synthesized using 10 transmit/ 


R, of 
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Figure 18. FFT beamforming with V,=3 knots and 


(a) R=100 yards, (b) R=250 yards, and (c) R=500 yards. 


The initial target angle, @=45°. The array is synthesized 
using 10 transmit/receive cycles. 
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Figure 18 shows the results of FFT beamforming with 


V,=3 knots and (a) R=100 yards, (b) R=250 yards, and (c) 





R=500 yards. The initial target angle, @=45° in all three 
cases and the array is synthesized using 10 trans-— 


mit/receive cycles. We again see that phase weighting the 








output electrical signals using the binomial approximation 





of received acoustic signal phase is almost exact. Phase 








weighting using the stop-and-hop approximation method for 


estimating received acoustic signal phase results in very 








low array gain in all cases. 


The errors introduced by the binomial approximation 





method seen in Figures 14, 15, and 16 are not large enough 


to significantly affect the steering and focusing of the 








beam pattern for typical vehicl speeds and target ranges. 








When the array is placed within 1 yard of the target and 


the vehicle is allowed to travel at 7 knots, however, the 





binomial approximation errors becom vident. Figure 19 








shows the results of FFT beamforming with R=1 yard, 


V,=7 knots, and @=45° for an array synthesized with 10 





transmit/receive cycles. In this case, the errors intro- 
duced in the estimation of propagation ranges with the 


binomial approximation method are a significant percentage 





of the total range. 
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Figure 19. FFT beamforming with R=1 yard, V,=7 knots, and 





@=45° for an array synthesized with 10 transmit/ 


receive cycles. 
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D. SUMMARY 


The propagation ranges and received acoustic signal 


phase can be estimated using both the stop-and-hop method 


and the binomial approximation method. 


duce some error into the 





and received acoustic signal phase. 





Both methods intro- 





calculation of propagation ranges 





indicate that the error introduced 





mation is significantly lower than 


the stop-and-hop approximation for 


speeds and target ranges. 


Simulation results 


by the binomial approxi- 





the error introduced by 





typical LMRS vehicle 





Estimates of received phase from 


both the stop-and-hop approximation and the binomial ap- 


proximation methods were 


electrical signals in an 





errors introduced by the 


used to phase weight the output 


FFT beamforming algorithm. 


The 








stop-and-hop approximation method 





significantly degrade th 








ing algorithm. 


method for received acoustic signal 


used to accurately phase weight the 


Conversely, 


performance of the FFT beamform- 


the binomial 


approximation 








nals at typical LMRS vehicl 


thesis conclusions are presented in 
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speeds 


phase estimation can be 





output electrical sig- 


and target ranges. The 


the next chapter. 
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V. SUMMARY AND CONCLUSIONS 


This chapter provides a summary of the propagation 
range and received acoustic signal phase estimation meth- 
ods. Based on the simulations presented in the last chap- 
ter, the accuracy of the stop-and-hop approximation and the 


binomial approximation methods is analyzed. Finally, fur- 





ther research in SAS acoustic signal processing is recom-— 


mended. 


A. SUMMARY 


Accurate estimation of acoustic signal propagation 








range and received phase is critical to SAS system perform— 





ance. To determine the phase of the acoustic signal inci- 





dent upon the receiver at all locations along the synthetic 


array, the time-varying range from the transmitter to the 








target and from the target to the receiver must be deter- 
mined. The estimated phase of the received acoustic signal 
can then be used to steer and focus the beam pattern of the 


synthetic aperture array. 


The SAS concept is currently being applied to the 








problem of water-borne mine classification and imaging in 





the Navy’s Long-Term Mine Reconnaissance System (LMRS). 





This requires high resolution imaging at long ranges. 


Under these conditions, small errors in the estimates for 





propagation range and subsequent phase estimates can sig- 








nificantly degrade the SAS system’s ability to steer and 


focus the array beam pattern. 


This thesis analyzed the performance of the stop-and- 


hop approximation and the binomial approximation models for 
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determining propagation range and received phase of acous-— 


tic signals for a SAS system imaging a stationary target. 





The stop-and-hop method, presented in [6], initially as- 





sumes that the vehicle is stationary during transmission, 
propagation, and reception of the acoustic signal. A sim- 


ple correction is applied to account for receiver motion. 





The binomial approximation method, presented in [5], devel- 








ops equations for the range from transmitter to target and 
from target to receiver for the general problem of bistatic 


scattering. Approximate solutions for propagation range 





equations are obtained by performing a binomial expansion 


of the propagation range equations. The derivations of 





this technique place some restrictions on the simulation 


parameters, but at the same tim nsure minimal estimation 








error. 


The stop-and-hop approximation method and the binomial 





approximation method were compared using simulations that 





exactly determine the acoustic signal propagation range and 








received phase for a SAS system imaging a stationary tar- 


get. This exact method, presented in [9], for the first 





time allows us to investigate the accuracy of common ap- 


proximations and their effect on SAS system performance. 


This thesis also presented an overview of the model 





used to simulate the propagation of acoustic signals for a 


SAS system. In the simulations presented, the ocean is 





treated as a linear, time-variant, space-variant filter. 





The complex frequency response of the ocean filter is the 


primary focus of the simulations presented. 
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B. CONCLUSIONS 
Simulations presented in this thesis are based on 


typical LMRS vehicle speeds and target ranges. Performance 





of the different approximation methods was compared for a 


Single transmit/receive cycle and then the estimates of 








received acoustic signal phase were used to compare the 
performance of each method applied to an FFT beamforming 


algorithm. 


The stop-and-hop approximation introduced a maximum of 


-0.5% received phase error with the vehicle traveling at 





the LMRS maximum specified speed and at the maximum target 





range. The error in propagation range and received phase 
estimation increased as target range increased, resulting 


in a maximum of -2.5% received phase estimation error with 











the transmit/receive platform at the minimum LMRS speed and 














the target at 25 yards. The estimation errors introduced 


by the stop-and-hop approximation method are small, but 








significant. This is best demonstrated when the stop-and- 











hop approximation method is used to determine the phase 





weights in an FFT beamforming algorithm applied to a SAS 


array. The error introduced by the stop-and-hop approxima- 








tion method significantly degraded the performance of the 








FFT beamforming algorithm. It was also evident that the 





accuracy of the stop-and-hop approximation method decreased 





as the target aspect angle increased. This suggests that 
the stop-and-hop approximation method is poorly suited to 
modern SAS techniques that combine multi-aspect data to 


improve image resolution. 


The binomial approximation method estimated acoustic 


signal propagation range and received phase per trans-— 
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mit/receive cycle with much less error than the stop-and- 





hop approximation method. At the maximum LMRS speed and 








target range, the received phase estimation error was ap- 





proximately 100 times less than the error introduced by the 


stop-and-hop approximation. Similarly, at low platform 








speeds and close target ranges the estimation error intro- 








duced by the binomial approximation method was 1000 times 


less than the error introduced by the stop-and-hop method. 








In fact, the errors introduced by the binomial approxima- 





tion were insignificant when the estimates of received 


acoustic signal phase were used to phase weight the output 





electrical signals of the array in an FFT beamforming algo- 





rithm. Only at very close target ranges (less than 1 yard) 








and high platform speeds did the errors introduced by the 





binomial approximation significantly affect the performance 








of the SAS system. In contrast to the stop-and-hop ap- 





proximation method, the accuracy of the binomial approxima-— 





tion method improved as the target aspect angle increased. 


The propagation range and received acoustic signal 





phase can be determined exactly for the SAS system using 


the techniques presented in [9]. These exact solutions do 





not add significant complexity to the SAS signal-processing 


problem while eliminating errors introduced by approximate 





methods. These exact solutions place no restrictions on 


target range and vehicle velocity, allowing for maximum 





flexibility as SAS systems evolve. 
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Cc. RECOMMENDATIONS FOR FUTURE RESEARCH 
This thesis considered improvements to SAS system 


beamforming using the exact and binomial approximation 





methods for determining propagation range and received 
acoustic signal phase versus using the stop-and-hop ap- 


proximation. The methods presented in [9] also provide 





exact solutions for the time-varying angles of incidence 
and scatter at the discrete point scatterer and the time- 
varying angles of incidence at the receiver. This addi- 


tional information should be incorporated into models of 








the scattering functions for mines and mine-like objects 





since scattering functions are not only functions of fre- 











quency, but also functions of the angles of incidence and 


scatter. Scattering functions should be included in future 








simulations in order to determine their impact on SAS sys- 








tem performance. Using LMRS operational parameters, typi- 
cal oceanographic data, and real world target simulations, 


and end-to-end SAS simulation should be developed. The 





simulations would be critical to the development of SAS 


operational guidelines, capabilities, and future improve- 





ment. 
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